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ABSTRACT 

MACHO-97-BLG-41 is a gravitational microlensing event produced by a lens composed of multiple masses 
detected by the first-generation lensing experiment. For the event, there exist two different interpretations of 
the lens from independent analyses based on two different data sets: one interpreted the event as produced 
by a circumbinary planetary system while the other explained the light curve with only a binary system by 
introducing orbital motion of the lens. According to the former interpretation, the lens would be not only the 
first planet detected via microlensing but also the first circumbinary planet ever detected. To resolve the issue 
using state-of-the-art analysis methods, we reanalyze the event based on the combined data used separately by 
the previous analyses. By considering various higher-order effects, we find that the orbiting binary-lens model 
provides a better fit than the circumbinary planet model with Ax^ ^ 166. The result signifies the importance 
of even and dense coverage of lensing light curves in the interpretation of events. 
Subject headings: gravitational lensing: micro - planetary systems - binaries: general 



1. INTRODUCTION 

The last two decades have witnessed tremendous progress 
in gravitational microlensing experiments. On the observa- 
tional side, improvement in both hardware and software has 
contributed to the great increase of the detection rate of lens- 
ing events from tens of events per year at the early stage of the 
experiments to thousands per year at the current stage. In ad- 
dition, photometry based on difference imaging substantially 
improved the quality of photometry. 

Along with the observational progress, there also has been 
advance on the analysis side. A good example is the analysis 
of light curves of lensing events produced by multiple masses. 
The light curve of a single-lens event is described by a sim- 
ple analytic equation with a small number of parameters and 
the lensing magnification varies smoothly with respect to the 
lensing parameters. As a result, observed light curves can be 
easily modeled by a simple minimization method. How- 
ever, when events are produced by multiple masses, modeling 
light curves becomes very complex not only because of the 
increased number of lensing parameters but also because of 
the non-linear variation of lensing magnification with respect 
to the parameters. The non-linearity of lensing magnifications 
is caused by the formation of caustics which denote positions 
on the source plane where the lensing magnification of a point 
source diverges. Caustics cause difficulties in lens modeling 
in two ways. First, they make it difficult to use a simple lin- 
earized minimization method in modeling light curves be- 
cause of the complexity of the parameter space caused by the 
singularity. Second, magnification computations for source 
positions on a caustic are numerically intensive. As a result, 
modeling a multiple lens event was a daunting task at the early 
stage of lensing experiments. However, with the introduction 
of efficient non-linear modeling methods such as the Markov 
Chain Monte Carlo (MCMC) algorithm, combined with ad- 
vances in computer technology such as computer clusters or 
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graphic processing units, precise and fast modeling became 
possible and now it is routin e to model light curves in real time 
as lensing events progress ("Dong et al.''2006'; 'Cassan"2008i 
Kains et al. 2009; Bennett 2010; Rvu et al. 2010; Bozza et al\ 
Furthermore, current modeling take into account vari- 
ous subtle higher-order effects. 

In this paper, we reanalyze the lensing event MACHO-97- 
BLG-41 that is a multiple-lens event detected by the first- 
generation lensi ng experiment Ma ssive Compact Halo Ob- 
jects (MACHO: :Alcock et anil993h . For the event, there ex- 
ist two different interpretations. Base d mainly on the data 
obtained from the MACHO experiment. iBennett et"an (Il999h 
interpreted the event as produced by a circumbinary planetary 
system whe re a planet was orbi ting a stellar binary. On the 
other hand, Alb row et al.l (2000), from independent analysis 
based on a different data set obtaine d by the Probing L ensing 
Anomalies NETwork (PLANET: lAlbrow et al .11 1998b group, 
arrived at a different interpretation that the light curve could 
be explained without the introduction of a planet but rather 
by considering the orbital motion of the binary lens. Accord- 
ing to the interpretation of Bennett et al. (1999), the lowest- 
mass component of the triple-mass lens would be not only 
the first planet detected via microlensing but also the first cir- 
cumbinary planet ever detected. Despite the importance of 
the event, the issue of its correct interpretation remains unre- 
solved. Therefore, we revisit MACHO-97-BLG-41, applying 
state-of-the-art analysis methods to the combined data used 
separately by B ennett et al. (1999) and Albrow et aL(2000,) . 

2. DATA 

The data used for our analysis come broadly from two 
streams. The first stream comes from MACHO survey obser- 
vations plus the Global Microlensing Alert Network (GMAN: 
Alcocketal. 1997) and Microlensing Planet Search (MPS: 
iRhie et al.iJ999i) follow-up observations. We refer to this data 
set as the "MACHO data". The other stream comes from ob- 
servations conducted by the PLANET group. We denote this 
data set as the "PLANET data". In Table [T] we present tele- 
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Fig. 1. — Light curve of MACHO-97-BLG-41 based on the combined MACHO plus PLANET data sets. Also presented is the best-fit model curve from our 
analysis. The two insets in the upper panel show the enlarged view of the two caustic-involved features at HJD' ~ 619 and 654. The lower three panels show the 
residual from the orbiting binary, triple, and static binary lens models. 

solution. During this normalization process, 3cr outliers from 
the solution are removed to minimize their effect on model- 



TABLE 1 
Data sets 



MACHO data 




PLANET data 




observatory 


number 


observatory 


number 


MS0 1.3m(R) 


711 


SAAO 1.0 m(/) 


97 


MS0 1.3m(B) 


772 


SAAO 1.0 m(V) 


14 


MSO 1.9 m(R) 


16 


ESO/Dutch 0.9 m (/) 


58 


Wise 1.0 m(R) 


17 


ESO/Dutch 0.9 m (V) 


18 


CTIO 0.9 m (R) 


35 


Canopus 1.0 m (/) 


95 






Canopus 1.0 m (V) 


14 






CTIO 0.9 m (/) 


49 






Perth 0.6 m (/) 


26 



Note. — MSO: Mount Stromlo Observatory; CTIO: Cerro Tololo Inter- 
American Observatory, SAAO: South Africa Astronomy Astronomical Ob- 
servatory. ESO: European Southem Observatory. 

scopes, passbands, and the number of data of the individual 
data sets. Note that Bennett et al. (1999) analyzed only the 
MACHO data set while , Albrow et al.. (.2000) conducted their 
analysis based only on the PLANET data set. 

In order to use data sets obtained from different observa- 
tories, we rescale error bars. For this, we first readjust er- 
ror bars so that the cumulative distribution of x^/dof ordered 
by magnifications matches to a standard cumulative distribu- 
tion of Gaussian erro rs by introducing a quadratic error term 
dBachelet e t alJl20T2l) . Then, error bars are rescaled so that 
X^/dof becomes unity, where is derived from the best-fit 



ing. 

In Figure [T] we present the light curve based on the com- 
bined MACHO + PLANET data. The light curve is character- 
ized by two sepai-ate peaks at HJD' = HJD-2450000 ~ 619 
and 654. We note th at the models of lBennett et al.1 (1 1 999 ) and 
lAlbrow et alj (120001) are consistent in the interpretation of the 
overall light curve including the peak at HJD' ~ 654. How- 
ever, the two models differ in the interpret ation of the peak at 
HJD' - 619. While Bennett et alJ m9§ ) explained the first 
peak by introduci ng an additional len s component of a cir- 
cumbinary planet, lAlbrow et al.l (120001) described the peak by 
considering the orbital motion of the binary lens. 

3. ANALYSIS 

3.L Standard Binary-lens Model 

To describe the observed light curve of the event, we be- 
gin with a standard binary-lens model. Basic description of a 
binary-lens event requires seven lensing parameters. The first 
three of these parameters describe the geometry of the lens- 
source approach including the time of the closest lens-source 
approach, fo, the lens-source separation (normalized by the 
angular Einstein radius of the lens, ^e) at that moment, mq, 
and the time scale for the source to cross the Einstein radius, 
fE (Einstein time scale). Another three parameters describe 
the binary nature of the lens including the projected separa- 
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TABLE 2 
Best- FIT Parameters 



parameter 


aLdllLldliJ. 


binary lens 


Ul UILtUcU clllaA 


■ 

triple lens 


y\ 1 1'ci 1 1 n "Y 
UtU alluA 


X^/dof 


2355.8/1920 


1915.1/1920 


1912.4/1920 


2086.6/1920 


2077.3/1920 


to (HJD') 


653.519±0.004 


653.426±0.006 


653.426±0.007 


653.371 ±0.008 


653.388±0.008 


MO (10-2) 


8.123±0.024 


-7.346±0.079 


7.303±0.074 


7.407±0.047 


7.284±0.049 


(days) 


20.26±0.016 


24. 11 ±0.229 


23.95±0.213 


24.82±0.155 


24.50±0.150 


^1 


0.494±0.0003 


0.479±0.002 


0.481±0.002 


0.485±0.002 


0.480±0.002 


?! 


0.481 ±0.003 


0.346±0.005 


0.342±0.005 


0.318±0.005 


0.326±0.005 


a 


1.211±0.003 


-1.176±0.003 


1.179±0.003 


1.152±0.003 


1.156±0.005 


i'2 








1.900±0.007 


1.875±0.008 


12 (10-3) 








5.149±0.160 


4.644±0.180 










1.872±0.003 


1.858±0.008 


P* (10-3) 


9.69±0.08 


7.49±0.12 


7.91±0.11 


7.36±0.09 


7.35±0.01 








-1.59±0.38 




0.21±0.15 








0.29±0.10 




0.34±0.08 


di/rff (yr-') 




-0.72±0.050 


-0.67±0.05 






<io/A (yr') 




-0.95±0.037 


-0.30±0.16 







tion s (in units of 9e) and the mass ratio q between the binary- 
lens components, and the source trajectory angle with respect 
to the binary axis a. The two peaks of the light curve of 
MACHO-97-BLG-41 are likely to be features involved with 
caustic crossings or approa ches during wh i ch finite-source 
effect becomes important ( DominikI 119951: iGaudi & Gouldl 
Il999t IGaudi & PettersI 120021) . To account for this effect, an 
additional parameter, the normalized source radius = 6** / 9e 
is needed, where 9^ is the angular source radius. In our stan- 
dard binary-lens modeling, we additionally consider the limb- 
darkening variation of the source star surface brightness by in- 
troducing linear-limb-darkening coefficients. Fx. The surface 
brightness profile is modeled as S\ oc 1 - F^Cl - 3cos(/)/2), 
where A denotes the observed passband and (p is the angle be- 
tween the normal to the surface of the source and the line of 
sight toward the source (Alb row et al. 1999). Based on the 
source type (subgiant) determined by the spectroscopic ob- 
servation conducted bv'L ennon et alj ( 119971) . we adopt coeffi- 
cients from Claret (2000). The adopted values are Tb = 0.793, 
Tv = 0.666, Tr = 0.575, and T, = 0.479 for data sets acquired 
with a standard filter system. For the MSO 1.3m data, which 
used a non-standard filter system, we adopt {rB + Tv)/2 for 
the B-band data and (Tr + Ti)/! for the 7?-band data. 

Although th e tw o sets of solution s pres ented by 
iBennettetall (119991) and l Albrow et al.] (12000 1) already 
exist, we separately search for solutions in the vast parameter 
space encompassing wide ranges of binary separations and 
mass ratios in order to check the possible existence of other 
solutions. From this, we find a unique solution with s ^ 0.49 
and q ~ 0.48. See Table [T| for the complete solution of the 
model. We note that the solution is b asically consistent w ith 
the results of lBennett etan(ll999h andlAlbrow et all (120001) in 
the interpretation of the main part of the light curve including 
the peak at HJD' ^ 654. At the bottom panel of Figure [T] we 
present the residual of the standard binary model. It is found 
that there exist some significant residuals for the standard 
model. This is also consistent with the previous analyses that 
a basic binary model is not adequate to precisely describe the 
Ught curve. 

3.2. Higher-order Effects 

The existence of residuals in the fit of the standard binary 
lens model suggests the need for considering higher-order ef- 
fects. We consider the following effects. 

First, we consider the effect of the motion of an observer 



caused by the orbital motion of the Earth around the Sun. This 
"parallax" effect causes the source trajectory to deviate from 
rectilinear, resulting in long-term deviations in lensing light 
curves (Gould 1992). The event MACHO 97-BLG-41 lasted 
~ 100 days, which is an important portion of the Earth's or- 
bital period, i.e. 1 year, and thus the parallax effect can be 
important. Considering the parallax effect requires two addi- 
tional parameters tte.n and tte.b, that are the two components 
of the lens parallax vector tte projected onto the sky along 
the north and east equatorial coordinates, respectively. The 
magnitude of the parallax vector corresponds to the relative 
lens-source parallax, TTiei = AU(Dl' - Z); '), scaled to the Ein- 
stein radius of the lens, i.e. tte = 7r,ei/6'E (lGouldll2004l) . 

Second, the orbital motion of a binary lens can a lso cause 
the s ource trajectory to deviate from rectilinear (iDominikI 
Ll99_8|). The orbital motion causes further deviations in lens- 
ing light curves by deforming the caustic over the course of 
the event. The "lens orbital" effect can be important for long 
time-scale events produced by close binary-lens events for 
which the event duration comprises an important portion of 
the orbital period of the lens system (Shin et al. 2013). To first 
order approximation, the lens orbital motion is described by 
two parameters, ds/dt and da/dt, that represent the change 
rates of the normalized bin ary separation and the source tra- 
jectory angle, respectively ("Albrow et al.ll2000l) . 

Third, we also check the possible existence of a third body 
in the lens system. Introducing an additional lens component 
requires three additional lensing parameters including the nor- 
malized projected separation, S2, and the mass ratio, q2, be- 
tween the primary and the third body, and the position angle 
of the third body with respect to the line connecting the pri- 
mary and secondary of the lens, V^. 

3.3. Result 

We test models considering various combinations of the 
higher-order effects. In Table |2] we list the goodness of the 
fits and the best-fit parameters for the individual tested mod- 
els. From the comparison of the models, we find the following 
results. 

First, we confirm the result of Albrow et al.l (1200 0) that the 
consideration of the orbital effect substantially improves the 
fit. We find that the improvement is Ax^ ~ 441 compared to 
the static binary model. When we additionally consider the 
parallax effect, the improvement of the fit, A^- ^ 2.7, is very 
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Fig. 2. — Geometry of the lens system for the best-fit binary (upper panels) and tiiple (lower panels) lens solutions. In each panel, the closed cuspy figures 
represent caustics and the curve with an arrow is the source trajectory. The filled circles is the locations of the lens components and the dotted circle is the Einstein 
ring centered at the barycenter of the lens. Right panels show the enlarged view of the coiTesponding shaded regions in the left panels. For the binary model, 
caustics and lens positions vary in time due to the orbital motion. We present two sets of caustics at HJD' = 619 and 654. The coordinates are co-rotating with 
the binary axis so that the binary axis aligns with the abscissa. All lengths are scaled by the Einstein radius Sg. 



meager. This implies that among the two effects, the orbital 
motion of the lens is the dominant higher-order effect in ex- 
plaining the residual from the standard model. In Figure [T] 
we present the model light curve on the top of the observed 
light curve and the residual of the model. In Figure|2] we also 
present the geometry of the lens system. 

Second, we find that the existence of a third body does 
not provide a fuUy acceptable fit. Our best-fit solution of 
three-body lens mo deling is consistent with the solution of 
iBennett et al.l ([1999') in the sense that the third body is a cir- 
cumbinary planet with a small mass ratio. See Table |2]for the 
best-fit parameters and Figure |2]for the lens system geometry. 
With the introduction of a planetary third body, the fit does 
improve from the standard binary model with Ax^ ^ 270. 
The additional consideration of the parallax effect further im- 
proves the fit with Ax^ ^ 9, but it is still substantially poorer 
than the orbiting binary-lens model with Ax^ ~ 166. From 
the comparison of the residual (see Figure [T), it is found that 
the triple lens model cannot precisely describe the light curve 
during and before the first peak, 595 < HJD' < 625. 

The key PLANET data that exclude the triple lens are from 
SAAO. Their smooth "parabolic" decline signals a caustic 
exit along the axis of cusp. This is compatible with the tri- 
angular caustic generated by the close binary, but not with the 
quadrilateral caustic induced by the putative planet. In partic- 



ular, the near-symmetric shape of the quadrilateral caustic im- 
plies that a cusp-axis trajectory would have a near-symmetric 
excess flux in the approach to the first peak as after its exit, 
which is not seen in the pre-peak MSO data. 

4. CONCLUSION 

We have conducted a reanalysis of the event MACHO-97- 
BLG-41 for which there exist two different interpretations. 
From the analysis considering various higher-order effects 
based on the combined data sets used separately by the previ- 
ous analyses, we find that the dominant effect for the deviation 
from the standard binary-lens model is the orbital motion of 
the binary lens. 

The result signifies the importance of even and dense cov- 
erage of lensing light curves for correct interpretation of grav- 
itational lenses. For MACHO-97-BLG-41, the difference be- 
tween the two previous interpretations partially stems from 
the poor coverage of the first peak that is important in the in- 
terpretation. Although a strategy based on survey plus follow- 
up observations can densely cover anomalies occurring at 
an expected time (e.g., peak of a high-magnification event) 
or long-lasting anomalies, it would be difficult to densely 
cover short-lasting anomalies arising abruptly at an unex- 
pected moment. Since the first-generation lensing experi- 
ments, there has been great progress in survey experiments. 
The cadence of survey observations has increased from ~1 - 
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2 per day to sever al dozens a day for the current lens ing ex - 
periments (OGLE :|Udals kil(l2003h. MOA : Bond et al. ( 12001b . 
ISumi et all (f2003h . Wise: IShvartzvald & Maoz. (,2012)). Fur- 
thermore, a new survey based on a network of multiple tele- 
scopes (KMTNet: Korea Microlensing Telescope Network) 
equipped with large format cameras is planned to achieve a 
cadence of more than 100 per day. With improved coverage. 



the characterization of microlenses by future surveys will be 
more accurate. 
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